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NOMENCLATURE 


C dimensionless parameter = (pp./p t p e ) 

C'jlf value of dimensionless mass transfer coefficient at zero mass 
injection rate 

Cjj value of dimensionless heat transfer coefficient at zero mass 
injection rate 

C'p species averaged specific heat, e.g. [=] BTU/lbm-deg R 
Cpi specific heat for pure species i 

c, constant relating thermal and mass diffusivity, (c.f. eq.(26) 
Ref.(9) 

0 reference binary diffusion coefficient(c.f. eq.(l9) Ref.(9)) 
e-g- [=] ft 2 / sec 

T thermal diffusivityyfor species i 

multi-component diffusion factor for species i (eq.(l9) Ref.(9)) 

static gas enthalpj', e.g. [=] BTU/lbm 

static enthalpy for species i 

species averaged gas enthalpy 

enthalpy of surface material 

enthalpy of surface material removed by melting or erosion 
total enthalpy (static plus kinetic energy) 
diffusive mass flux of element k, e.g. [=] lbm//t 2 -sec 
mass fraction of element k 
t molecular weight of gas mixture 
Mi molecular weight of species i 
m c mass removal rate of surface material by direct 
vaporization processes, e.g. [=] lbm//f 2 -sec 
>h g mass injection gas rate from any in-depth pyrolysis process 
T i mass removal rate from surface by melting, etc. 

total pressure, e.g. [=] atm. 
r species averaged Prandtl number 


v 


q c diffusive heat flux, e.g. [=] BTU//f 2 -sec 
Qcand heal conducted internally to the heated surface 
q r radiative heat flux to the body surface 
R universal gas constant. 

r local streamline radius measured normal to the body symmetry axis, 

e -6- [=] ft 

r 0 local body radius 

Sc species averaged Schmidt number 
T local fluid temperature, e.g. [=] deg. R 
u boudary layer tangential velocity, e.g. [=] ft/sec 
U e boundary layer edge tangential velocity 

v boundary layer nomal velocity component 

x, mole fraction of species i 

Greek Symbols 

a*., mass fraction of element k in species i 

/3s nose radius parameter (c.f. eq(23), Goulard Ref.(ll)) 

X gas phase thermal conductivity, e.g. [=] BTU/ft-sec deg R 
p gas phase viscositj', e.g. [=] lbm/ft-sec 
4>k elemental reaction source term, e.g. [=] lbm//f 3 -sec 
p gas phase density, e.g. [=] lbm//t 3 

Subscripts and Superscripts 

e boundary layer edge value 
FC fully catalytic state 
i species i 
k elemental species k 
N C non catalytic state 
v local radius exponent, 

)' = 1, for axi symmetric 
i' = 0, for planar 

tii body surface or ’’wall” value 
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ABSTRACT 


Transpiration cooling of hypersonic aircraft and planetary entry vehicles has 
often been proposed to help reduce surface heating and maintain operating 
temperatures within the limits of available materials technology. For exam- 
ple, proposed aircraft which can take off horizontally and directly achieve 
orbital velocity (Single Stage to Orbit, (SSTO)) will experience very high 
heating rates under flight conditions of high enthalpy and high stagnation 
pressure. The dominant surface heating mechanism for these situations 
will be conduction and convection with relatively low levels of radiative 
heating. Transpiration cooling by mass addition can be quite effective in 
blocking such dissipative heat fluxes. This paper presents a quantitative 
study of the effect of transpiration cooling on a specific SSTO maximum 
heating case. The analysis is performed for the stagnation region of the 
vehicle using both similarity theory and a numerical, non-similar reactive 
boundary layer solution program, BLIMPK. An outline of both methods of 
analysis is presented, and comparisons are made between the determined 
heat fluxes. Additional results are presented for the effect of mass addi- 
tion of pure nitrogen, helium and argon on heat flux and boundary layer 
species profiles. The effect of transpiration on the heat flux is examined for 
surfaces of different chemical reactivity. It is shown how the non-similar 
numerical, BLIMPK, results can be used to develop a heat transfer coeffi- 
cient correlation in the presence of mass addition from the analytical heat 
flux expressions obtained from similarity theory. Suggestions are made for 
future research and experimental study. 
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I. INTRODUCTION AND BACKGROUND 


Transpiration cooling by injection of fluids into hypersonic boundary layers 
has long been recognized as an effective means of mitigating the effects of 
conductive and convective heat fluxes on high speed re-entry vehicles. (Ref. 
1,2,3 ) It is also a primary mechanism in blocking dissipative heat fluxes 
in ablative heat shield systems. As such, this active means of thermal pro- 
tection will be most efficient in situations where the dominant surface heat 
flux is dissipative rather than radiative. Usually transpiration cooling would 
then be applicable for use on nose cones and flight (and control) surface 
leading edges of relatively small curvature typical of high L/D re-useable, 
re-entry vehicles. The flight regime of these applications will often pass 
through regions of maximum total heating at quite high stagnation pres- 
sures. As a result the dominant heating mechanism will be conduction and 
convection. This paper presents an analysis of the transport phenomena 
associated with transpiration into hypersonic boundary layers on generic 
blunt bodies of revolution (sphere-cones). Results will be presented for a 
specific trajectory maximum heating point of a hypersonic, Single Stage 
to Orbit (SSTO) vehicle. In the future, such aircraft are envisioned to 
take off horizontally from a runway and continuously thrust directly into 
orbit. As described by M. E. Tauber (Ref. 4) these flights will encounter 
maximum stagnation heating during the ascent phase, and this will be al- 
most entirely dissipative in nature. Optimum design and determination of 
thermal protection materials requires knowledge of surface heat fluxes and 
temperatures. Transpiration cooling can dramatically affect these quanti- 
ties, and this study will provide some quantitative estimates of the effects 
of operating parameters on surface heating for an SSTO ascent. 

An additional phenomenon which can significantly affect heating in hy- 
personic flight is the reactivity of vehicle surface material. When atomic 
species from dissociated boundary layer gases reach the relatively cool sur- 
face of the vehicle, they can recombine and release their formation energy 
directly to the surface. The amount of energy released will depend on the 
inherent surface reactivity (or catalycity) and the rate of supply of atoms 
to the surface. Transpiration can most definitely impinge on the rate of 
atom supply, and this interaction will be examined. 
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II. SCOPE OF WORK AND OBJECTIVES 

The objective of this sudy is to develop an approximate methodology ca- 
pable of calculating the degree of transpiration cooling possible at various 
fluid injection rates near the stagnation region of high Mach Number exter- 
nal flows such as occur on SSTO type vehicles. Analytical solutions derived 
from similarity assumptions as well as solutions obtained from a numeri- 
cal non-similar boundary layer computer code (BLIMPK) will be applied. 
Information obtained here will be useful in determining candidate thermal 
protection materials. Additional detailed information on the effect of mass 
injection rate and surface catalycity on boundary layer atomic species pro- 
files will be examined and related to the effect of the transpiration cooling 
mechanism. Finally, surface heat flux expressions derived from similarity 
theory will be compared with the numerical, non-similar results. This will 
provide a limited check on the validity of the simplified method. 

III. DEVELOPMENT OF APPROXIMATE ANALYSIS 

The mathematical analysis of mass injection into reactive, hypersonic shock 
layers requires a fundamental approach, necessitated by the alteration of 
the surface boundary conditions. Basically the addition of mass at the 
boundary introduces an inhomogeneity which changes the shape of the 
boundary layer velocity, enthalpy and species concentration profiles. Even 
without this complication, a complete two or three dimensional solution 
of the thin shock or boundary layer governing equations requires complex 
numerical methods. It is, however, still important to ascertain if certain 
approximate methods will suffice to provide engineering estimates of surface 
heating. 

Attempts at analytical description of mass injection into boundary lay- 
ers began at least as early as 1951 with W. E. Stewart (Ref. 5) of M.I.T. He 
presented a similarity solution to the case of heat and mass transfer from 
the surface of a heated (or cooled) fixed temperature flat plate immersed 
in a boundary layer flow with mass addition at the surface. A constant 
property generalized solution was presented. This work will be alluded to 
later in this paper. In 1961 and somewhat later P. A. Libby and others 
(Ref. 6,7 ) obtained numerically augmented similarity solutions to the case 
of stagnation region transpiration cooling of hypersonic blunt bodies at 
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high mass injection rates. In 1959 P. M. Chung of NASA Ames Research 
Center (Ref. 6 ) supplied a partial similarity solution to stagnation region 
mass injection in the presence of finite rate surface reactions. J. T. Howe 
(Ref. 2) has given a synopsis of correlations, semi-analytical solutions, and 
numerical solutions for high Mach Number stagnation region mass addi- 
tion problems (c.a. 1964). All of these approaches have in common the 
presumption of reactive thermo-chemical equilibrium among the boundary 
layer dissociation reactions as well as equilibrium at the shock or boundary 
layer edge. Many features of these previous attempts will be shared by the 
approximate methods and results presented in this paper. 

III.l THE GENERAL TWO-DIMENSIONAL, MULTICOMPONENT 
REACTIVE BOUNDARY LAYER PROBLEM WITH MASS 
ADDITION AT THE SURFACE 

There are a very large number of authors to solutions of this general prob- 
lem, and no attempt is made to refer to them in this study. Problably 
the most lucid and still perhaps the most accurate solution of the two- di- 
mensional axisymmetric problem is the integral matrix method numerical 
solution of E. P. Bartlett and R. M. Kendall (Ref. 9 , 1965). This method, 
embodied in the computer code BLIMPK, has been adapted to the spe- 
cific transpiration cooling problems examined here, and reacting boundary 
layer results from this code are given below. A brief description of the code 
and solution characteristics is also provided. At this point, however, the 
boundary layer governing equations of Bartlett and Kendall will be used as 
a basis set for the level of description desired in this study. These equations 
will also serve as a starting point for the more simplified analysis presented 
later. 


III.l. 1 GOVERNING EQUATIONS 

In the continuum, hypersonic limit under the assumption that the shock and 
boundary layers are thin compared to the body dimension, the following 
conservation equations are applicable for a two-dimensional, axisymmetric 
blunt body. The coordinate system (s,y) is curvilinear and fixed to the body 
surface with origin at the stagnation point; (s) is the streamline variable, 
and (y) is the body normal variable. 
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Total Continuity: 


1 d(pur t/ ) + djpv) _ Q 


r o 9s 


dy 


Conservation of Species (i): 


dK k dK k d 
P U ~K~ + P V ~T~ ~ ~a~ 0 *) + 4>k 


ds 


dy dy 


where, 


3k = 


pDp_ 2 

Pi M 


dZ k 


dy 


+ (Z k — A*) 


dlnui 


dy 


+ Y1 a ki^i 


T dlnT 


dy 


Streamwise Momentum Conservation: 


dv 


du d ( du\ dP 


P% ds + pX dy dy \ r dy 
Body Normal Momentum Conservation: 

pu 2 dP 


ds 


dy' 


Conservation of Total Energy: 


dH T dH T d 
' , ” - 5r + '”"s7 = Ty {q - + q ' ] ' 


where, 


d(ir) dT pDpi n 


where, 


(1) 


( 2 ) 


(3) 


(4) 


(5) 


( 6 ) 


(7) 


«• “ | - (*>■ + S) f + c ‘* T t + C ' RT "*) If (e| 
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These equations are to be supplemented by mechanical and caloric equa- 
tions of state as well as the appropriate boundary conditions. Specifically, 
Hj is the total enthalpy (static plus kinetic) 


o 


H t = h + J 

(9) 

h = y / K i h, 

(10) 

hi = f C Pi dT + h° 

(11) 

Cp = Y ] KiCpi 
| 

(12) 

K k = Ta k ,K, 

i 

(13) 

Zk = Y2 a kiZi 

( 14 ) 

Zi = MiXi/FiPi 

(15) 

mi = Y2 x i F i 
% 

(16) 

= 2- F . 
i 

(17) 


(IS) 


(19) 

The simplest boundary conditions in effect here 

are the following; 

u(s, 0) = 0 

(20) 

p(s,0)t>(s,0) = p w v u , 

(21) 

h(s, 0) = h w 

(22) 

K k {s,0) — K kw . 

(23) 

Equations (22) and (23) could alternatively be replaced by wall diffussive 

flux balances. 
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III. 1.2 TRANSFORMATION OF VARIABLES. 


The above partial differential governing equations are, in general, non-linear 
and parabolic in form. A hypersonic blunt body has a sub-sonic flow region 
in the vicinity of the stagnation point, and an iterative problem must be 
solved if the proper boundary conditions are to be imposed on the above 
equations. In thin shock layer and boundary layer problems, this is usually 
accomplished in an approximate manner. Inviscid shock layer governing 
equations can be solved to provide an approximation of the body surface 
pressure distribution. With this as a fixed boundary condition the bound- 
ary layer system is amenable to various numerical and semi-analytic so- 
lution techniques. In taking this approach here, the BLIMPK computer 
simulations with surface mass addition represent the parabolic numerical 
solutions, while similarity analysis will be used to obtain simpler analytic 
results. In either method, the first step is to reduce the order of term 
by term and variable coupling effects through the use of variable trans- 
formations. The transformation of variables will group dominant effects 
into specific terms. In axisymmetric boundary layer theory the Levy-Lees 
transformations are the most useful. These are; 


f = fp.U.^ds 

(24) 

* <v. /» . 
p v 

(25) 

II 

fc- 

-3 

(26) 


is the stream function. Omitting algebraic details for clarity, (see Ref. 9 ) 
the boundary layer equations transform over to the following; 


Continuity: 


pv = -■ 


Pt l'ePe r o 


V'-s 


01 

' +(2 ‘>4j 


- A/5f-' 

rf ^ Os 


(271 


7 



Species (i) Conservation: 


fK'u + I — + 4>k 


A f 


dK 




dln£ k dln{ ) 

£ Momentum: 

"" + vrr - % m - «• = 2 (/'J£ - r$) 


V' = — + P' 


at ( >/5? l A. 


1 ? Momentum: 


P'= l -AM(ff 




Energy Conservation: 


w +(-*.+ *r - 1 


(28) 


(29) 


(30) 


(31) 


(32) 


The individual diffusive fluxes annotated in the above equations take the 
form; 

( C \ rz, 


Jk = 




(33) 


P4 = In (P2 T c ‘) 


and, 


■?o = “ 




oe 

S 


/7"f ; t 2 + + n' 

rv 


(34) 


(35) 


w ' = i 


;V - (t> + — )r + c,jwv; + (h-h + c ( rt H3 ) p; 
V Pl/'2/ 


(36) 


8 


The boundary conditions for the transformed system can take on a variety 
of different acceptable values, but the following axe most typical. 



f edge jj 1 edge 

(37) 

Cge = 


(38) 

Ht\ edge 

= Hx(actual) e j ge 

(39) 


Hip Ue = 0 

(40) 

^k\edge = 


(41) 


a-;u. = 0 

(42) 
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(43) 


#r(f>0) = Htw 

(44) 


K k ({,0) = K kw 

(45) 


Htw a»d Kuw can also be indirectly determined by the following surface 
mass ans energy balances. Whether these are fixed values or subject to the 
constraint equations depends on the actual physical nature of the surface 
boundary. 

™ a h g + m c h c - £ m T ,hi - (p w v w )h w - q cw + q TU , - q corui = 0, (46) 

i 

nigKkg 4- ™ c K kc - V m r ,K k i - (p w v w ) K kw - j kw = 0. (47) 

i 

111.2 NON-SIMILAR SOLUTIONS AND BLIMPK 

The terms on the right hand sides of the transformed governing equations 
and asscociated terms in the dissipative fluxes are known as non-similar 
terms and can be significant in situations where either the surface or bound- 
ary layer edge conditions have large variations along the streamwise di- 
rection. This can happen with abrupt changes in body geometry, sudden 
changes in surface conditions such as abrupt changes in mass addition rates 
or surface reactivity, or rapid changes in edge conditions due to external 
flow effects. 

Non-similar effects can only be evaluated by solving the complete set 
of governing equations. The current version of BLIMPK has been adapted 
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to numerically solve for the general case of mass addition and surface fi- 
nite rate chemistry. BLIMPK uses an integral matrix procedure augmented 
with Newton-Raphson methods to numerically solve the non-similar bound- 
ary layer equations. In this study, solutions are obtained for the injection 
of pure nitrogen, helium and argon in the stagnation region of an SSTO 
vehicle. The surface of the vehicle is considered to be adiabatic and at a 
temperature corresponding to a surface energy balance where all incoming 
conductive and chemical energy is fully re-radiated to the vacuum. Fluid 
mass injection rates can be specified independently at any streamwise body 
station, and surface temperatures, species mass fractions and fluxes can be 
calculated. Plots of these quantities are shown and described in Section IV. 

III.3 SIMILARITY ANALYSIS AND SOLUTIONS FOR 
TRANSPIRATION WITH SURFACE REACTION 

A great many situations in hypersonics either do not exibit non-similar 
behavior or these effects can be considered to be second order. In the 
vacinity of the stagnation point of blunt bodies, the boundary layer will 
be essentially similar in the streamwise direction. The calculation of sur- 
face temperatures and dissipative fluxes in this region is important in an 
engineering sense, since the stagnation region usually experiences the most 
severe combined heating in most flight situations. It will be useful to obtain 
approximate analytical expressions for these quantities and compare them 
with the more general non-similar results from BLIMPK. 

To obtain a simplified solution it is assumed that all terms containing 
explicit dependency on the streamwise variable, f , are small and can be 
dropped. The more general boundary layer model outlined in Section III. 1 
includes the effects of both finite rate, homogeneous gas dissociation reac- 
tions as well as kinetic surface catalyzed re-combination reactions. Non- 
equilibrium kinetics effects are very non-linear and must be handled nu- 
merically. A great many cases of hypersonic blunt body flow will have high 
enough shock layer pressures to allow the dissociation reactions to reach 
thermochemical equilibrium by the time the flow’ reaches a boundary layer 
edge. This would be especially true in the subsonic region. Starting from an 
equilibrium edge, driving forces for displacement from chemical equilibrium 
are almost non-existent within the boundary layer, except near the surface 
where catalysis can be taking place. Under these conditions the reaction 
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source terms can also be eliminated from the governing equations. This 
assumption is almost universal in the similarity treatment of mass addition 
problems. Including these caveats, a set of similar boundary layer equa- 
tions can now be derived. For this situation where the total shock layer 
thickness is considered very thin, the body normal momentum equation 
will be dropped, (i.e. P' — 0). The transformed similar boundary layer 
equations now become, 


Total Continuity: 


PtV'H'K x 

pv = ~^r f 


Conservation of k-th Elemental Species: 

Sc 




Streamwise, £ , Momentum: 


Energy: 


Pr 

C 


"r\ f ~ (£)] h t + h t = c (i - 1) v\rr 


Boundary Conditions will now take the form; 


(48) 


(49) 


(50) 


(51) 


/(0) = fun /'( 0) = 0, f'(oo) = 1 (52) 

Ht(0) = H w , Ht(oo) = H' t (53) 

A'fc(O) = A7, K k (oo) = Kl (54) 

The right hand sides of the momentum and energy equations still contain 
inhomogeneities associated with coupling of conduction, viscous dissapation 
and compression effects. As Lees (Ref.10) has pointed out, these terms are 
second order also for most hypersonic flows, since Pr == 1 and f' 2 — p t /p ~ 0 
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when T w T t . These assumptions are usually referred to as the ’’Cold 
Wall” approximation. With this premise then, a closed form solution is 
possible. W. E. Stewart (Ref.5 ) has previously obtained this solution in 
the form, 





(55) 


satisfying the generic governing equation, 


where, 


ft 

c 


-(§)]* 


+ A" = 0, 


A(0,ft, —fw) = 0 

A(oo,ft, —f w ) = 1 


and, 



1.0 

Sc 

Pr 


Formally the solution is, 


A(»), ft, —fw) 


/>*?(- r§( 

'/-(£)')<*>) 

1 dr) 

c ‘*r (- /: g 

(/-(§)')*»! 

) di ) 


(56) 

(57) 

(58) 

(59) 

(60) 


with 

A'( 0 ,n,-/ w ) = j />P (-£§(/ - (I)') *) *»}" («» 

Using this kind of analysis it is desired to obtain the total heat flux to the 
surface of the body in terms of the solution variables. This can be written 

Z 1 S 

Qw = 9u. + ^2jk^k(Pw) (6-) 

k 

where, 

X = p.v.ch (ft - &:) < 63 > 
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C M = C. | r (^)’ / "| (Sc”)" ft'IO, A*-,-/.) (64) 


and 


qu, = f> e U t C H (HZ-H‘ T ) 

r -r li£.\ 

Cm {pt w ) Ki'(0,Sc\-f u 


(65) 

( 66 ) 


'Sc w \ H t '( 0, Pr w , —fw) 

,) 

After performing the appropriate algebra, the following explicit expression 
for can be obtained. 


<?: = (Ca/Ci ) n- T + v (C ^ C&) 

Jk ^rxn 


7r<=: 


-r (ft - ft 8 ) + (ft - ft) 


where, 


4'“' = 1 - 

x ram x 


Wi = 


e.v.Cu 7 r*r 

f)\U filD 

(Hj — Hj ) 


*r 


nr 

(67) 

( 68 ) 

(69) 

(70) 

pJ7.Cft(JI*-jy*<.Ftf)) (/1) 

Two mathematical limits of tins equation are of interest. These refer to 
the cases of either fully reactive (catalytic) or completely non-reactive wall 
reactions. Thus, for a fully catalytic wall k%-+oc> , for all k, 

QUFC) = (CVCJ) + £<£WCS) (Ki - ft 8 ) HT(FC) (72) 

k 

and for a non-catalytic wall, AiJ"— +0 , 

Q-JNC) = A W«J(.VO + , ..’"ft- - ft ^ (ft - A-;) Wr(JVC) (73) 

( 711 U, ~ A( a/ ) fc C w 


r (Hf(FC)-Hi.) 

nr = 

Ql = 


(Hjt(FC) - Hj) 

(-<£) 
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where, 


Pu,V w 

^ " PeU c C° M 

(74) 

A fa = Ch/Cr 

(75) 

■A/a i = Csi/Clf. 

(76) 


Goulard (Ref.ll) has obtained the additional special limiting case for a cold 
wall, stagnation point value for C° H • 

lim P'U'Ch = y/2p ef i e /3 N (. Pr w )~ * fo.47 (Pr“) 1/3 | . (77) 

tJv — 0 \ Pepe / t J 

Equations (67), (72) and (73) have a well recognized behavior as a function 
of «?*,. A/# and A I'm are exponentially decreasing in rii ^ (see Section V), and 
this leads to the general behavior as all previously published plots for mass 
addition with a fully catalytic wall. It shows the classic form of the blowing 
correction function for ablation and transpiration. However, the intermedi- 
ate and non-catalytic cases exibit a peak at lower values of mass injection. 
This non-calssical behavior is due to addition of the injected gas enthalpy, 
which will be greater than the dissipative energy blocked by injection at 
low injection rates. This will occur, most probably, only in those situations 
where the injected gas equilibrates rapidly with the boundary layer. If 
there is a significant cold gas layer heat transfer resistance, then the added 
gas enthalpy will be less, and the maximum will be less or non-existent. 
Ail added feature of this maximal behavior displayed at low catalycity is 
the singular behavior of Equation (73). The maximum must occur when 
rh*,— vV A /. For the heat flux to remain finite at this point, , in 

other words, the mass diffusive heat flux must become zero. This will be 
the balance point between heat flux due to mass diffusion and injected mass 
flux. This is consistent with intuition. 

IV. COMPARATIVE NEAR-STAGNATION SOLUTIONS FOR 
DIFFERENT GASEOUS INJECTANTS 

In addition to formal theoretical development, an objective of this study is 
the determination of quantitative transpiration cooling effects from stag- 
nation region gas injection for a specific hypersonic flight case. Using the 
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non-similar numerical simulation, BLIMPK, a series of calculations were 
performed to compare the relative effectiveness of different molecular weight 
gases for stagnation cooling of an SSTO during maximum heating. As men- 
tioned above, in hypersonic cruise, the maximum heating point will occur 
during ascent to orbit. Flight conditions typical of this are shown in Table 
1. (Ref. 4 ) BLIMPK solutions were obtained for a range of stagnation 
point injection mass fluxes using nitrogen, helium and argon as coolant 
gases. An additional parameter in these calculations was the surface re- 
activity. Resultant heat fluxes correspond to a fully adiabatic vehicle wall 
with all energy ultimately being re-radiated to the vacuum. The effect 
of surface reactivity can be seen from the results for an equilibrium edge 
calculation for this small radius (1.5 inch) blunt nosed cone configuration. 
A plot of surface temperature and streamline coordinate is shown in Fig- 
ure 1 . Curves are given for a relatively non-catalytic silica tile material 
coated with borosilicate glass, and for a case where the surface reactivity 
was increased 56-fold at a distance 2.5 inches from the stagnation point. The 
additional chemical energy available from a catalytic surface is significant, 
and in this SSTO case the stagnation heating is severe; even for a non- 
catalytic surface. Stagnation wall temperatures as high as 5000 deg F are 
not feasible for state-of-the-art high temperature materials. One hope of 
transpiration cooling is to bring these high surface temperatures down to 
within the safe operating range of existing materials. 

Table 1 

Flight Conditions for Maximum SSTO Ascent Heating 

Freestream Conditions: Stagnation Conditions: 

Pressure: 1.1091 x 10 -3 atm. 0.924 atm. 

Density: 9.0475 x 10 -5 lbm/ft3 
Velocity: 26,380 ft/sec 
Altitude: 155,000 ft. 

Enthalpy: 13,899 BTU/lbm. 
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As mentioned above, a comparison study of different molecular weight 
gases as injectants was performed using BLIMPK. Figure 2 shows a plot 
of stagnation heat flux vs nitrogen injection rate for the sharp nose cone 
portion of the SSTO at maximum heating. This plot is dimensional to give 
an intuitive feel of the magnitude of required injected gas. For this flight 
condition, the freestream mass flux, > is approximately 1 lbm//f 2 - 

sec; therefore, the injection rates are effectively dimensionless relative to 
the freestream. With reference to Figure 2 it is seen there exists the same 
form of functional behavior as the similarity solutions embodied in quations 
(67), (72) and (73). It is encouraging then, that for constant mass injection 
rates into laminar boundary layers, BLIMPK could be used to provide a 
semi-analytical correlation of the heat and mass transfer coefficient expres- 
sions and the parameters contained therein. Figures 3 and 4 show plots 
of stagnation transpiration cooled heat fluxes for helium and argon injec- 
tion respectively. The relative influence of molecular weight is dramatic 
and in favor of lower weight species such as helium and hydrogen. This is 
further shown in Figure 5 where the fully catalytic results are shown on a 
common plot. Also shown here are the interpreted results of Chung (Ref. 
8), which occur as straight lines for his fully catalytic and non-catalytic 
cases. The linear behavior is due to maintenance of a constant value of the 
function, ( 4*“ xn = 0 and 1.0) for variable injection rate. This, of course, is 
not occurring for the plots derived in this study, but the trend is the same. 
Energy release to the surface by catalytic recombination is important, and 
the effect of mass injection on this is discernable from the calculations in 
this study. The derived functionalities of total heat flux vs mass injection 
rate as well as the BLIMPK calculations show a diminishing return for the 
non-catalytic surface at. higher and higher injection rates. An examination 
of the competing phenomena at the surface will shed some light on why 
this occurs. Figures 6 and 7 are BLIMPK results for species mass fraction 
profiles across the boundary layer at the stagnation point, respectively for 
a fully catalytic and non-catalytic surface. These results are for a zero 
mass injection rate and show that the boundary layer (with equilibrium 
edge) is essentially frozen at equilibrium except for the elemental nitrogen 
for the fully catalytic surface. It is also noted that the surface catalyzed 
gas phase dissociation reaction kinetic model considered here, oxygen total 
surface recombination is much less active than nitrogen. This is consistent 
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with the higher activation energy for nitrogen dissociation. Figures 8 and 9 
show the same BLIMPK results for the very high rate of nitrogen injection 
of 0.30 lbm//f 2 -sec. The large mass addition of nitrogen at the surface 
has flattened out the species profiles at the surface, and ail profiles appear 
as almost identical. In the non-catalytic case, even though the elemental 
nitrogen no longer reacts at the surface, it is replaced at the surface by the 
injection of molecular nitrogen. As a result the overall effect is nearly the 
same as for a fully catalytic surface at these same high injection rates. 

V. DISCUSSION OF RESULTS AND EVALUATION OF SIMILARITY 

THEORY 

The functional results from similarity theory and BLIMPK discussed above 
are interesting and useful for estimating the order of magnitude of transpi- 
ration cooling needed to bring surface operating temperatures to within the 
operating limits of materials used for the thermal protection of hypersonic 
vehicles. Such estimates, at this point, can be considered only as order- 
of-magnitude in quality. Limitations which must be imposed include the 
existence of turbulence and shock layer interaction with injected gas. Injec- 
tion of high mass fluxes in the stagnation region is known to trip turbulence 
at very low Reynolds numbers, resulting in fully developed turbulence just 
immediately off the stagnation point. If this occurs, the above estimates 
of surface heat flux as a function of mass rate of injection will be too opti- 
mistic. There is very little knowledge of the transition levels for this type 
of boundary layer flow. When injection rates get even higher than those 
resulting in turbulence, interaction with the bow shock will occur. The 
degree of shock displacement, or even the character or existence of a well 
defined shock will be in question. An entirely more general analysis would 
then be required to calculate heat fluxes. 

Recalling that Figure 2 and equations (67), (72) and (73) indicate a 
good correspondence between similarity theory and the more general non- 
similar BLIMPK results near the stagnation point, it is possible to use the 
numerical data from the non-similar model to calculate specific correlational 
forms for the flux expressions obtained in similarity theory. If, as is plausible 
for most cases, it is assumed that C'm/Cjj=1.0 , then equation (67) contains 
only two adjustable parameters or functions. These are, 

A h = Ch/Cjj ("*>) 
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and, 


= Cm/CIj. ( 79 ) 

Many previous studies have, e.g. Ref.(12), have shown that the typical 
functionality for M can be represented as. 


M = 


Km: 


(e Kth » - 1)’ 


( 80 ) 


where, k can be considered an adjustable constant, k is usually of the 
order of unity, Ref.(l2). This functional form is responsible for the classic 
exponentially decaying behavior for the surface conductive heat flux as a 
function of injected mass flux. This empirical form can then be used in place 
of the BLIMPK calculations where it can be assumed that the boundary 
layer maintains a self-similar or constant pattern form over the flight body. 

VI. CONCLUSIONS AND RECOMMENDATIONS. 


In flight situations where -the freestream enthalpy and stagnation pressures 
are relatively high, the dominant surface heat flux will be due to conduc- 
tive/ convective (dissipative) mechanisms. Transpiration cooling appears 
to offer an effective means of blocking this energy flux in the same manner 
as occurs in ablating systems. It is observed in this study that the beneficial 
effect of a non-catalytic surface becomes negligible at high mass injection 
rates. At these high mass rates, net surface heat transfer will be quite low, 
however. The absolute value of required injection rate for nearly complete 
heat blockeage is uncertain due to a lack of detailed knowledge of transpi- 
ration on boundary layer turbulence and shock interference. It does appear 
that, for stagnation region heat transfer, similarity theory is an adequate 
description of mass addition in the presence of finite rate surface reactions. 

Future studies should include definitive experiments to ascertain the 
injection rates required to trip turbulence for a variety of body shapes 
and flight conditions. Numerical computational studies should also be per- 
formed to examine the interaction between injected gas and the bow shock 
at very high injection rates. The threshold for this phenomena should also 
be investigated experimentally in an arc-jet environment. All of these ex- 
periments should range over enough conditions to be comparable to the 
above calculations and any future results. 
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Figure 1 


SSTO Surface Temperature Profiles 

showing effect of surface reactivity 
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Figure 2 


SSTO Surface Heat Flux vs. Wall Mass Injection 
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Figure 5 


SSTO Surface Heat Flux vs. Wall Mass Injection 

Comparison of Fully Catalytic Results 
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Figure 6 


SSTO Stagnation Boundary Layer Species Profiles 

Full - Catalytic N2 Injection Mdot = 0.0 Lb/ft2— sec 

□ = N 
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Figure 7 


SSTO Stagnation Boundary Layer Species Profiles 

Non-Catalytic N2 Injection Mdot = 0.00 Lb/ft2-sec 
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A«K0 
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Figure 8 


SSTO Stagnation Boundary Layer Species Profiles 
Full - Catalytic N2 Injection Mdot = 0.30 Lb/ft2-sec 

□ = N 
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A-NO 
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Figure 9 


SSTO Stagnation Boundary Layer Species Profiles 

Non — Catalytic N2 Injection Mdot = 0.30 Lb/ft2— sec 

□ = N 
0-0 
A-NO 
+ -02 
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